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ABSTRACT 


Six different methods for averaging historical salinity 
values were studied for their applicability in the determin- 
ation of depth corrections due to sound speed variations. 

The averaging techniques range in complexity from the use of 
aesingle value of salinity throughout the entire water column 
to the use of a historical salinity average profile corrected 
for the surface salinity. Results show that historically 
determined salinity can be used in the computation of the 
@epen Corrections without exceeding the accuracy limits. Two 
mmere Values of salinity were found to be sufficient to cover 
the West Coast of the United States: 31 °%., applicable north 


mameooN, and 33 /, for use south of that latitude. 
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NERO OUT TON 


fee lie ECHOSOUNDER 

Despite the development of laser- and photobathymetry, 
the echosounder, or fathometer, is still the most commonly 
used method for the determination of depth in hydrography. 
The transducer emits a sound pulse which is reflected by 
the bottom and received back at the transducer. The time 
of travel of this sound pulse is divided by two and this 
value multiplied by the assumed speed of sound in sea-water,- 
meomeiving the depth according to the expression z = vt 
(distance = speed x time). This transformation 1s made 
either mechanically or electronically within the fathometer 
Meseit, and the result displayed is the nominal or fatho- 
meter depth beneath the transducer. This depth is not equal 
tO the true depth, since the assumed sound speed generally 
does not equal the true speed of sound in sea-water, which 
varies throughout the water column. One example of an 
observed sound speed profile is presented in Figure l. 
For determining the true depth, one must know the sound speed 
profile throughout the water column and apply it according to 


the following expression [Greenberg and Sweers, 1972]: 


Peenough sound speed 1s the correct term, most hydro- 
graphic literature uses sound velocity instead. 
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some echosounders have the possibility of being set to 
different sound speeds, in order to represent as closely as 
possible the average speed of sound in the water column. In 
feet cases, however, they are constructed with a calibrated 
speed of sound assumed in the device, usually 1463.04 or 
1500 m/s (800 or 820 fathoms per second). The echosounders 
Msed by the National Ocean Survey (NOS) are calibrated for 
an assumed speed of sound of 1463.04 m/s which, according to 
Umbach [1976] is a value reasonably close to the sound speed 
mmmrest Waters surveyed. Corrections are then applied to 
miemaepths obtained by the fathometer in order to compensate 
for the use of a sound speed other than the one that actually 


affects the travel of the sound ounces 


peeeeesCCURACY REQUIREMENTS 
The International Hydrographic Bureau (IHB) established 
in 1968 the following requirements for accuracy of depth 


determination [International Hydrographic Bureau, 1968]: 





; micnough these are depth corrections, the term usually 


applied in hydrographic literature is sound velocity correc- 
irons . 





Depth Range pecuGacy skequUlLEemene 


wetl fathoms (0 - 20 m) eters 
me55 fathoms (20 - 100 m) SO) it 
deeper than 55 fathoms Poot eines depth 


The above standards represent the maximum allowable error in 
miemm@easurement of depth due to all sources, including the 
Variations of sound speed within the water column. When con- 
Mieerineg Crrors due to sound speed alone, no errors can be 
Saiameeto Or exceed 0.257 of the depth, which means that the 
mean sound speed must be known to within +4 m/s [Umbach, 1976}. 
The standards mentioned above apply only to conventional echo- 
sounding systems, in which the sound pulse beam is vertical. 
Mener systems, applying oblique acoustic paths, require higher 


accuracy standards. 


fee rED OF SOUND IN THE SEA 

Speed of sound in a fluid depends on the density (9), and 
Compressibility (k) of the fluid and the ratio between the 
specific heats of the fluid at constant pressure and at con- 
Stant volume (Y), according to the following expression 


[Sverdrup, Johnson and Fleming, 1942; Bowditch, 1972]: 


eee 
Ok 
As these parameters depend on temperature, salinity and 


pressure (which depends on depth), it is more convenient 
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to measure these properties instead, and compute the sound 
speed from them. 

Different expressions have been developed empirically 
by Kuwahara [1939], Del Grosso [1950], Wilson [1960a, b], 
Leroy [1969] and, more recently, Clay and Medwin [1977]. 
Mackenzie [1960] compared several of these formulas, which 
showed differences between them of less than 3 m/s, or 0.2% 
meer o depth [Urick, 1975]. Speeds used for depth correc- 
tions in NOS are calculated using Wilson's equation [1960b], 
which is described in Appendix A. This equation is the most 
commonly accepted formula due to its large range of applica- 
peeeey and its relative accuracy (0.30 m/s standard deviation) 


[Wilson, 1960b, 1962]. 


eee RADITIONAL METHODS FOR DETERMINING DEPTH CORRECTIONS 

Several different procedures can be used to determine 
Mmmemaepth corrections, or correctors, to be applied to echo- 
soundings. These corrections can either be found directly 
or by knowing the sound speed profile. The sound speed can 
either be measured directly (direct comparison systems) or 
computed after the determination of some characteristics of 
the water column, as stated above (oceanographic and limno- 
logic determination). Some of the different procedures are 
described below: 


l. Direct Comparison Systems 


The depth recorded by the echosounder is compared 


to the actual depth indicated by the leadline or the depth 


eo 





at which the bar is lowered. Because these are comparative 
methods, not only are sound speed variations included in the 
results but also instrument errors [Umbach, 1976]. Both 
methods are limited to shallow depths and good weather con- 
Gieeions. 

a. Leadline Comparison (Vertical Cast) 

The leadline is lowered to the bottom at several 
depths over a range which should be at least as great as that 
likely to be encountered in the survey area [The Hydrographer 
of the Navy, 1969]. Several readings are taken from the 
fathometer at each depth, and all are plotted individually 
meorder to eliminate random errors. In areas where launch 
hydrography joins or overlaps ship hydrography, it is of 
particular importance to use this method on both ships and 
launches in order to resolve possible discrepancies in the 
soundings [Umbach, 1976]. This method is limited to calm 
seas, shallow depths, and areas of hard bottom for accurate 
readings. 

b. Bar Check 

The depth of a bar indicated by the echosounder 
1s compared to the known depth at which the bar is lowered 
below the transducer. The procedure is used at several depths 
and the readings are taken both when the bar is lowered and 
when it is raised. Since the bar is handled over the side, 
this method is usually limited to small vessels and shallow 


depths, although it is sometimes performed mechanically on 
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Mmeeee vessels. The depth limitation still holds, however, 
due to difficulties in maintaining vertical lowering lines 
faeein keeping the bar within the beam of the transducer. 
mieeeprocedure should only be used in calm seas with light 
Mees and no strong current. 
fee Direct Measurement of Sound Speed 

Sound speed can be measured directly by lowering a 
sound velocimeter to various depths and recording the speeds 
indicated. Although accuracies obtained by this method 
can be of the order of a few tenths of a meter per second 
[Ulonska, 1972], sound velocimeters are not generally used 
because they are very expensive instruments [Calder, 1975]. 

Seeicasurements of Salinity, Temperature and Depth 

These parameters can be determined by direct meas- 
urement either with Nansen casts (temperature and salinity 
at specific depths) or with continuous sampling instruments 
which measure salinity (or conductivity), temperature and 
depth (STD or CTD). The latter instruments can also be of 
an expendable type. Sound speed at various depths is then 
computed using an equation relating salinity, temperature 
and depth (pressure) to sound speed, such as Wilson's equa- 
tion [1960b], or they can be extracted from tables, wherein 
sound speed is tabulated for different values of temperature, 
pressure and salinity [Heck and Service, 1924; Matthews, 


ieee, 19359; Kuwahara, 1939]. 


sy 





‘weedistouical Tables or Atlases 

Sound speed or salinity and temperature data have 
been tabulated for various areas and seasons. Matthews' 
memes | 1939] provide corrections to be applied to the depths 
obtained with fathometers calibrated to a sound speed of 1465 
or 1500 m/s. These depth corrections are given as a function 
Of nominal depth and are tabulated for 52 different oceanic 
areas identified by Matthews [Greenberg and Sweers, 1972]. 
feeommethod 1S not used by NOS but it is used in some coun- 
Beaes tor depths exceeding 200 m. For example, this is the 
ome thod which Australia uses for determining corrections 


ieweter depths of more than 200 m [Compton, 1980}. 


feet LAYER METHOD 

im Options 

In addition to the above techniques to establish the 

depth correction, yet another more sophisticated method is in 
widespread use--the "summation of layers" method [Umbach, 
1976]. This method is applicable when the sound speed is 
determined either by direct measurement or by computation from 
the depth, temperature and salinity parameters of the water 
column. This method can be applied in two ways, which are 


described in the "Hydroplot/Hydrolog Systems Manual"? 


ie il: 


[ Wallace, 





5 
The Hydroplot/Hydrolog system is an automated data aqui- 


Sition and processing system in use by NOS, respectively with 
and without a plotting capability. A description of this 
system 1S given in the Hydrographic Manual, Fourth Ed. [Umbach, 
m7). 
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a. Curve-fit Method 
The eee Solum tsmpanecitioned into layers of 
memeane thickness which are predetermined by the operator. 
Layer mid-depths are then computed and the sound speed at 
those depths interpolated either graphically or mathemati- 
Seely etrom a curve passed through the points at which the 
sound speed has been computed. 
Bee NO-cCUurve-f1t Method 
The water column is partitioned as above but the 
depths at which the observations have been made are consid- 
ered as the mid-depths, the boundaries between layers being 
defined as half-way between two successive observations of 
Mm@emcast. One example illustrating the no-curve-fit method 
is shown in Figure 2. 
ae computation of the Depth Corrections 
a. Layer Sound Speed 
For each layer a sound speed is determined for 
the corresponding mid-depth by any of the methods described 
above. This value of sound speed is assumed to apply through- 
omemene entire layer. 
mo. Layer Factor 
A factor is computed for each layer, based on the 
sound speed for which the instrument is calibrated, using the 
following expression [Umbach, 1976]: 
hea GC 
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Figure 2. Layer distribution for the no-curve-fit 


option of the layer method. 
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where A; is the actual sound speed for the i-th layer which 
has been determined for the corresponding mid-depth (see 
preceding paragraph), and C is the assumed sound speed of 
Eeewinstrument. 
mao Lavyer Correction 

The above factor is multiplied by the layer 
mime kness , Az. (ehewdatecrence between the depth at the 
boieeom and at the top of the layer) in order to obtain the 
momereetiOn to the depth (LC. ) due to the variation in sound 


Seed Of that layer: 


[iC = omeac tOKne x AZ. 
al 1 af 


Mmeemwater COlumn above the transducer is eliminated when 
determining the corrections for the upper layer. 
d. Bottom-of-layer Correction 
The layer corrections are added from the surface 
downward in order to obtain the total correction at the bottom 


of each layer: 


where Cc. is the correction at the bottom of the i-th layer. 
e. Fathometer Depth Correction Table 
The above corrections correspond to the true 
depth of the bottom of each layer. To obtain the fathometer 
depth to which these corrections will apply, one must sub- 


tract the correction (C.) from the true depth of the bottom 


aA 





of the corresponding layer Cae 


Cee 4 
1 


where Cen is the fathometer depth corresponding to the i-th 
1 


ayer. 
Pr CcOrreceLon Curve 

The values of the corrections are plotted against 
the fathometer depth. A curve is drawn to fit these depth- 
correction pairs. From this curve the depth correction appli- 
cable to any reading of the fathometer (within the range for 
which the curve has been computed) can be obtained. An example 
of such a curve, for which a transducer at zero depth has been 


considered, is shown in Figure 3. 


F, VARIATION OF THE PARAMETERS 

ies Variation of sound speed with depth, temperature 
Pgmecalinity depends on the initial conditions of these par- 
ameters. For a temperature of 14°C, salinity of 34 %o and at 
zero depth (sound speed under these conditions is 1502.94 m/s) 
an increase of one degree in temperature will increase sound 
m@eeed by about 4.5 m/s; an increase of 1 ¥ in salinity will 
increase the sound speed by about 1.3 m/s; an increase in 
depth of 100 m (about 10 atmospheres in pressure) will in- 
Beeasc SOund speed by about 1.7 m/s [Urick, 1975]. Latitude 
also has some influence, since pressure increases from the 


equator to the poles because of the increase of the gravita- 


tional attraction. This influence, however, is very small 


ae 
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and can be ignored. An example of the variation of sound 
speed with temperature, salinity and depth is shown in 
mpoure 1. 

Temperature is the parameter that has the largest influ- 
ence on sound speed. Also, it shows more variability with 
meme. Umbach [1976] states that, for the sound speed to be 
meeurate to within +4 m/s, temperature and salinity must be 
determined with an accuracy of +1°C and 1%», respectively. 
However, recent studies have been performed by the Testing 
Division, Office of Marine Technology (OMT), NOAA, regarding 
Mmiesintluence of the different parameters on sound speed. 
These studies show that if temperatures can be determined 
Beeer an accuracy of +.1°C, the salinity accuracy need only 


be within +3% to still meet the requisite accuracy standards 


femvins, 1976). 


Ge USE OF HISTORICAL INFORMATION 

Historical information has been used for some time in the 
computation of depth corrections for echosoundings, but usually 
only for depths in excess of 200 m [{Sherwood, 1974; Urick, 
mo). Yeager {1979] studied whether depth corrections of 
sufficient accuracy could be determined from historical salin- 
ity and temperature only. He chose an area which is repre- 
Semtative of the near shore region along the East Coast of 
the United States. The historical sound speed data there 
indicated that the temporal and spatial variability of tem- 


perature exceeded the acceptable limits for sounding 
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@oemrections, thus precluding the use of historical data alone 
for correction determination. In his study he also concluded 
mages in that area (a little more than one degree square with 
depths less than 200 m) salinity values were sufficiently 
stable to allow acceptable depth corrections to be computed 
meomenistorically derived salinity values combined with in 
facus temperature measurements. 

In Australia data have been collected over the contin- 
ental shelf since 1965 which show salinity to vary generally 
between 34.5 and 36 %o, changing only slightly with depth 
down to 300 m [Calder, 1975]. At present, the Australian 
Hydrographic Service exclusively uses a value of 35 Yo for 
the calculation of sound speed in depths less than 200 m 


feompton, 1980]. 


Bee OBJECTIVES 

ihe purpose of this thesis is to analyze the possibility 
of using historical salinity data averages for a large area 
with in situ temperature measurements to compute depth cor- 
rections and still remain within the required accuracy limits. 
The historical salinity data were treated using several dif- 
ferent techniques which increased in accuracy as more complex 
averaging schemes were used. 

Temperature-salinity observations generally have to be 
made at least once a month and be of sufficient number to 


cover the entire area sounded [Umbach, 1976]. The use of 
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historical salinity data and XBT (expendable bathythermo- 
graph) temperature information could save a significant 
amount of time and money if results remain within the 
accuracy standards. Also, since little time is needed for 
each XBT observation, more frequent temperature measurements 
could be made (both in time and space), thus increasing the 


accuracy of the corrections used throughout the survey. 
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Vereen OCe Uncen METHODOLOGY 


Some of the terms to be used in the following chapters 


are defined below: 


Historical data: data resultant from some averaging 
Scheme applied to previous observations in a large 
area. 


Inferred data: data obtained from both historical and 
observed salinities, used to represent the data for a 
Bance region. It has not been measured, and consists 
of a single value. 

Station: group of measurements taken at a particular 
Mice in a particular location. Each station consists 
of several observations taken over a given depth 
mocerval. 


Observation: measurement taken at a particular depth 
mone Station. 


matemcorrection: depth correction resultant from the 
application of the observed data. 


Estimated correction: depth correction for which the 
computation of some historical data were used. 


Error: difference between the estimated and true 
corrections. 


memeentage of error: division of the error by the 
depth over which the correction applies. 
A. STUDY AREA 
mae area selected for this study was the continental shelf 
region of the West Coast of the United States, extending off- 
eaore generally to depths of about 1500 m. The study area 
Samepe seen in Figure 4. This area was chosen on the basis 


of several factors: 
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Figure 4. Study area, showing the inshore limits of the 
Alaskan and Californian water masses. 
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1. It represents a fairly large area wherein a consider- 
Paeceamount of data were available. For the purpose of this 
meta a large area was desired to generalize the applica- 
Dility of the results. 

Eee Most hydrographic surveys are conducted over the 
Semeainental shelf region. 

3. The area was well suited for acquisition of histor- 
meeecdatd, as Will be seen later in this chapter. 

fee lhe extension to 1500 m depths guaranteed that all 
Shelf features, including submarine canyons, would be con- 
sidered. 

The study area was subdivided into three sub-areas: 
Pemenern, Southern and Strait of Juan De Fuca and Puget Sound. 
This latter area was considered separately because of its 
Megeacular salinity characteristics. The North-South division 
of the Pacific Ocean waters along latitude 40°N corresponds 
to the separation of the two different water masses consid- 


Soa in the historical salinity file. 


B. DATA SOURCES AND TREATMENT OF DATA 
Meets cOnmical Data 
Historical data for this study were obtained through 
the PROFIL subroutine of the Integrated Command ASW Predic- 
tion System (ICAPS). The ICAPS water mass history file was 
developed by the U.S. Naval Oceanographic Office (NAVOCEANO) 


orecOnecains historical information on the water masses of 
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the Northern Hemisphere and Indian Ocean [Department of the 
Navy, 1978]. Each ocean is divided into several areas and 
Subareas. In each of these subareas up to five different 
water masses may be included. Figure 5 shows the North 
Pacific Ocean, Pacific Area A divisions and the water masses 
present in Area A. The water mass files are also grouped 
seasonally: winter (January through March), spring (April 
through June), summer (July through September), fall (October 
mmeeiwteh December). 

Within the study area the ICAPS water mass file con- 
sists of only two water masses: Alaskan water (north of 40°N 
and Californian water (south of 40°N). These are shown in 
Figure 5(b). The inshore limits of these water masses are 
Shown in Figure 4. The Alaskan water mass was used for the 
Mear shore area off Washington, although this area is not 
mieluded in either of the two water masses. 

POmenGArSudata retrieval and processing, several dif- 
ferent program modules are available [Chace, 1979]. In this 
study only the PROFIL subroutine was used. One of the out- 
Pues Of PROFIL is historical salinity. It was used to obtain 
information for each area during all four seasons. The 
salinity profiles from these eight runs were punched onto 
cards and are shown in Table I. 

2. Observed Data 


The main source of observed data was the Fleet Numer- 


1cal Oceanography Center (FNOC) file entitled MOODS (Master 
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T2900 (°C) vT_(°Cc) 


Region Water Mass Name Hin Max Hin Max Poettion Freq. (2% 
Al TRANSL TION 7 1) l 35 
KUROSIIEO 13 19 2 65 
A2 FAST TRANSITION 10 16 1 100 
Aj NORPAC 7 32 1 43 
FAST TRANSITION 12 146 2 57 
AA CALI PORNLAN 5 11 1 60 
FAST TRANSITION il 19 2 40 
AS CALIFORNIAN 5 tl 1 ran 
A6 GULF OF CALIFORNIA 19 19 1 ran 
A? ALFUTLAN 0 7 1 8 
NORPAC ? 12 2 92 
Aa ALAS AN 0 ? I 44 
NORPAC ? 12 2 54 
Ad ALASV.AN 0 9 I 199 


AIO AL FUTIAN 0 4 I 190 
Figure 5(c) 
mere S(a) North Pacific Ocean locator chart, (b) Pacific 


area A, and (c) Pacific area A water masses, as defined 
miethe ICAPS history file. 
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TABLE I 


FROM ICAPS 


HISTORICAL SALINITY FOR THE WEST COAST 
OF THE UNITED STATES, 
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Oceanographic Observation Data Set). Some additional data 
were obtained from the Pacific Marine Center (PMC), NOAA, 
Seattle, Washington. The MOODS file contains data from sea 
water observations all over the world. Different types of 
observations are included: STD, CTD, Nansen cast, BT or XBT, 
and others. Each station contains the geographical location, 
date, time and a code for the ship/agency that made the ob- 
servation. Temperature, salinity, sound speed, oxygen con- 
Gentration or other properties are also provided depending 
on the type of observation. The depth interval between 
observations for different stations is not standardized. 

The time interval chosen for the study was the ten 
year period from 1964 to 1973. A more recent period was not 
chosen because a very large percentage of the recent data is 
from XBT observations, which does not provide salinity infor- 
mation. 

A direct extract from the MOODS file was not readily 
usable on the Naval Postgraduate School computer since the 
MOODS data are stored with a variable block size. A nine 
meack tape, 1600 bits per inch, constant block size, was then 
prepared at FNOC which contained temperature and salinity for 
Mmepens Standardized to 0, 30, 60, 90, 120, 150, 175, 200, 250 
and 300 m. Values in excess of 300 m were not given, even if 
they had existed in the original data. In the meantime, some 
data were punched onto cards, extracted from a printout ob- 


tained from the original data. Since this data set had a 
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feeeer protile definition, especially useful for the inter- 
mediate values between the surface and 30 m, it was decided 
momuse this punched data in substitution for the correspond- 
ies tations on the tape. 
A total of 3459 stations were analyzed, with the 

Mester obution shown in Table II. For each of these stations 
a depth correction was calculated based on the procedures 
mmerned im the previous chapter. Such corrections are con- 
sidered true depth corrections and are the values against 
which the estimated depth corrections (based on some mean 


historical salinity value) must be compared. 


cep enol 


DISTRIBUTION OF THE STATIONS BY AREA AND SEASON 


Winter Spring Summer Fall Total 
North 248 679 643 195 1765 
South 450 55 0 229 1953 Ze 


Puget Sound/ 
eenait of Ley 149 100 96 472 
Juan De Fuca ee 


ona | 3459 


Appendix B shows all the station locations for each 
area and season. Multiple occupations of a station may have 


been made which are not indicated by these figures. 
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C. DATA ANALYSIS 
ie oalinity 

From the ICAPS data file, mean salinity values were 
computed from the surface to a given depth, for each area 
and season. This mean was determined not only for the depth 
of each given value of historical salinity (Table I) but also 
for the mid-depths between two successive values. Table III 
contains the values obtained. These data were then used to 
establish a variety of historical salinity averages to be 
described below. 

Estimated depth corrections were then determined for 
each station using the observed temperature and an historical 
salinity value derived from one of four different methods: 

a. Gross Historical Salinity Average (Method 1) 

This method uses a single value of salinity 
throughout the entire water column. This salinity value 
was determined by averaging the historical salinity data 
from the surface to 200, 1000 or 3000 m. The choice of the 
depth depended on whether the deepest depth for the station 
Beeetess than 200, between 200 m and 3000 m, or greater than 
3000 m. 

b. Fine Historical Salinity Average (Method 2) 

This method was similar to the previous one 
Peeept that, instead of three possible values for salinity, 
ten possible values were considered: averages from the sur- 


face to 50, 100, 150, 200, 300, 500, 750, 1000, 1500 and 
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eeem.depth. As before, the depth greater than or equal to 
the deepest depth observed was chosen. For example, if the 
@epen were 250 m, the salinity average would be the one cor- 
responding to the 300 m water column. 

@eeeiltstorical salinity Curve (Method 3) 

The value from the historical salinity curve, 
which corresponded to a given depth, was obtained by interpo- 
itenon trom the existing values extracted from the ICAPS file. 

d. Adjusted Historical Salinity Curve (Method 4) 

Since the salinity value observed at the sea sur- 
face can sometimes deviate considerably from the historical 
surface value, an alternative to the historical salinity curve 
was considered. A linear adjustment was made to this curve 
so that the surface value would equal the observed salinity 
value, and the value at 30 m would equal the historical data 
(TCAPS) at that depth. Below 30 m the historical curve was 
used. 

Salinity values illustrating these first four methods 
memappiied to a particular station are shown in Figure 6. In 
light of the results obtained using these four methods, two 
new approaches were taken: 

e. Inferred Value (Method 5) 

A single value of salinity was used for each area, 
independent of season and depth. This value, which was inferred 
meome both the historical data and the actual observations, was 


judged to represent average conditions for the area. 
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meeure 6. Historical salinity derived from ICAPS. 


Ee 58 


Pe iHeukecamitmcmith Observed Suritace Salinity 
(Method 6) 


PiicomircenoCmiomswiitan tO the previous method, 
meeepe that the salinity which is applied to the first layer 
is the observed surface salinity. This method is also simi- 
lar to method 4, except it only uses one value of salinity 
throughout the column for depths other than the surface. 

meee COrrections 
Each of the six methods was applied as in a normal 

Survey. Estimated depth corrections were computed for each 
Meyereusing the no-curve-f1t method, i.e., the depth of the 
observed temperature was used as the mid-depth for the layer. 
iteemprocedures indicated in section I.E. were followed using 
steps a to d: 

(a) Estimated sound speed for each layer. 

(b) Computation of the corresponding factor. 

mee Layer correction. 


(d) Bottom of layer correction 


The correction curves for the observed temperature and salin- 
ity values and those based upon each of the first four methods 
are shown in Figure 7. The estimated correction at the bottom 
of each layer was compared to the "true" correction which was 


obtained using the observed salinity. 


Error, = estimated correction, - true correction. 
| il 


The differences were divided by the depth, in order to obtain 
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Figure 7. Correction curves for the first four methods. 
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miempercentage of error for that particular depth: 


ET ror. 
ae 


Meneentage of error, = : 
1 


where Z. momene deptheorethe bottom of the layer. If this 
value were greater than or equal to 0.25%, the method was 
mememaered Unsuitable at that depth, since the accuracy re- 
maprements could not be met. The values compared were the 
©Orrections at the true depth at the bottom of each layer. 
immereality, the values at the fathometer depth (true depth 
minus correction) should have been compared. This was not 
done since the fathometer depths would probably be different 
fOr each error. However, the resultant error from this approx- 
imation iS insignificant--of the order of a few millimeters 
(see Figure 8). 

MunEneNCOMmputations Lor the first four methods were 
executed sequentially in the same program, shown in Appendix 
C. This program was also used for the fifth and sixth methods. 
An output of this program for a particular station is shown in 
Figure 9. The meaning of the different headings is as follows: 

SH1: gross historical salinity average (method 1) 


SH2: fine historical salinity average (method 2) 


Mees depth of the observation 

ID : depth at the bottom of each layer 
T : observed temperature 

fae Observed salinity 
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fim historical salinity value (method 3) 

fips adjusted historical salinity curve (method 4) 
V : sound speed computed from Z, T and $ 

V1-4: sound speed computed using methods 1-4 


eee cCOrrection applied at depth ID computed 
from V 


CORI1-4: correction obtained using methods 1-4 


DCORI-4: error of methods 1-4 (difference between 
COR1-4 and COR) 


Pel - 4; BeGeceiuagcenOunecrnon Of methods 1-4 
(Ohi divaded by ID) 
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CPOs ONeOr ThE RESULTS 


Based on the foregoing analysis, four different possi- 
bilities have been considered for each of the six salinity 
averaging methods tested: 

(1) No error exceeds the limits at any given depth. 


(2) The error limit is exceeded only in the surface 
ayer. 


(3) The error limit is exceeded between the surface 
layer and the layer which contains the 30 m depth. 


Wry Excessive errors have been found for layers deeper 
than 30 m. 
Peeeeeoe OF HISTORICAL SALINITY DATA DERIVED FROM ICAPS 
(FIRST FOUR METHODS) 
ihe results obtained from using historical salinity data 
derived from the ICAPS water mass file are shown in Table IV. 
A discussion of these results, categorized by area, follows. 
1. Northern Area 
The large variation of surface salinity in this 
region caused the depth correction for the first layer to 
exceed the acceptable limits for the first three salinity 
averaging methods considered. Excessive errors were found 
for 31.6%, 28.0% and 24.4% of the stations examined, respec- 
tively. <A slight improvement is noted between methods 1 and 
3, indicating that some benefit is derived from the use of 


an historical mean salinity profile. However, when the 
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TABLE IV 


ERRORS RESULTANT FROM THE USE OF HISTORICAL SALINITY 
Dok eve D SEROMOICAPS (METHODS 1 TO 4). 
Results are shown for the surface layer and for water 
depths shallower or deeper than 30 m. Numbers indicate 
the number of stations which exceeded the accuracy 


standards. 


Percentages are indicated in parenthesis. 


ae ee Method 1 Method 2 Method 3 Method 4 
‘aa ee sor weto Sila <0) 330 4 Sale 30) «3002S SLL. Sc 30 >30 
Winter 248 93 38 3 78 24 2 73 2h 2 0 0 0 
(GSeecnmcls sly e515)  (9.7)> (0,8) (29,4) (8.5) (0.8) (0.0) (0.0) (0.0) 
Soring 679 301 100 44 278 151 31 235-1 3S 25 0 0 0 
eM ecoesy 4-9) et 2252) as. 6)) (54-6)) (19.9) .(3.7). (0.0) (0.0) (0.0) 
North saawer o83 15] So } 130 27 0 115 21 0 0 0 0 
Pao ctnCn mms) 20s (42) (On (17.9) 105.5) (0.0) (0.0) (0.0). (0-0) 
Fala Les Lv 3 Q 4 3 Q 7 3 0 0 0 0 
(Seletiesseo-0) mM.) (1.5) 000.0) (3.0) (1.5) (0.0) (0.0) (0.0) (0.0) 
fomieies) S57 tC f0imeeS. S35 80 CC for eo 0 
(Simoyeiiss sic. 7) 028.0) (ilo) (1.9) 024.4) (10/2) (1.5) (0.0) (0.0) (0.0) 
Winter 450 1 UV (} 0 0 0 Q) () 0 0 0 0 
(0.2.0) (Ol (0.0) 90010) (020) (0.0) 9 (0.0) (0.0) (0.0) €0.0) (.0) 
Spring 350 0 0 0 0 0 U 0) J 0 0 a Oi 
(GietOnG) "Ute  (0.0)) (OnO}a e(O.0) (0.0) “(O.0) (0.0) (0.0) (0.0) (0.0) 
South Sumer 225 3 0 0 1 Q 0 ( 0 0 0 0 0 
(1.3) (0.0) (0.0) (0.4) (0.0) (0.0) (0.0) (u.0) (0.0) (0.0) (0.6) (0.0) 
Fall 193 0 0 1 0 0 0 0 0 0 0 0 0 
COnuyee (Orgy. C0LS) (0.0) (0.0) (0.0) (0.0) (0.0) (6.0) (0.0) (0.0) (0.0) 
mae UC The 0. oO vO 0 0 ov O 
(0-3) (0.0) (0.1) (OoT)) 90.0)" (0.0) CH0)}s 40-0) (0.0) (0.0) (0.9) 
winter 127 blo 103 be) 102 Vo 2 35 yy 7y 0 1 1S 
Hoge el WC tHEU) AB SIn(7S.6) (7204). (os.d aoe. 2} (59.3) (0.0) (0.3) (21.8) 
Spring 149 131 129-125 126 125 121 125 123 L19 0 i 5 
(87.9) (86.6) (83.9) (84.0) (83.9) (81.2) (83.9) (82 9) (79.9) © (0.6) (0.7) (2-0) 
+ - 
Poor omer 100 v1 69 64 7U 68 o4 70 os a 0 2 8 
Pe be fuca (71.0) (69.0) (64.0) (70.0) (68.0) (04.0) (70.0) (03.0) 97.0) (Q.0) (2.0) (8.0) 
Fall yo 45 6d 53 S8 S6 54 47 Xe 34 v 4 + 
(o7.7) (60.7) (40.4) (00.4) (58.3) (50.3) ‘49.u) (3? S) (35.4) (0.0) (4.2) (4.2) 
Totul 472 377 365. S46 356 S45 331 328 Wr we |6hUMUlClU CO 
(79. eye croes sn (75.4) (73240 70.1) (CRRA 63.4) (60.6) (u.0) (1.7) (0.4) 
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historical salinity profile was adjusted to match the surface 
value (method 4), there were no stations which exceeded the 
acceptable limits. 

The depth error of the first layer was often so large 
that sometimes the offset created induced an unacceptably 
meee error in the second layer. This occurred for 14.9%, 
11.6% and 10.2% of the stations examined, respectively. In 
very few cases (2.7%, 1.9% and 1.5%) the excessive difference 
in surface salinities created an excessive error which extended 
to the third layer. 

There are two main reasons why this error is so large: 
1) the observed salinity extends to very low values, and 2) 
most of the data only are defined at the surface and 30 m 
isopleths which creates a very wide first layer (15 m) over 
which that value is applied. The errors which occur due to 
low surface salinities generally do not cause the accuracy 
limits to be exceeded except close to the surface, for depths 
rarely exceeding 30 m. The salinity may be very low at the 
feerace but it increases quite rapidly to nominal oceanic 
values in a few meters. 

The error due to low surface salinities is of par- 
ticular importance during the spring season (April-June) due 
to the large outflow of the rivers of Northern California and 
Oregon, especially the Columbia River. It is also noticeable 


in the winter season, but is considerably less important in 


summer. An example where this situation frequently occurred 
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was in the area of influence of the Columbia River, where 
on one occasion a salinity as low as 8.02°f0 was observed 
during the spring season. 
2. southern Area 

Considering the least sophisticated of the salinity 
averaging schemes (method 1), only five stations out of 1222 
(0.4%) showed any kind of deviation which exceeded the accept- 
able depth limits. Of these five, only one station, located 
adjacent to the Golden Gate Bridge in San Francisco, exceeded 
those limits in a layer other than the surface. For this 
weueren, Unacceptable depth errors were found at depths 
between 70 and 80 m due to an inversion in the salinity pro- 
file. Of the other four stations, only one exceeded the 
Famers limit when the second method was used. All these sta- 
tions were located in water more than 300 m deep, except one 
located about 8 miles south of Point Reyes, about half-way 
to the Farallon Islands. None of these observations exceeded 
the limits for the third and fourth methods. The general 
pattern resulted in a slight improvement in accuracy from the 
first to the third methods, and a significant improvement 
when the fourth method was applied. However, there were cases 
where, because of the large range in salinity, method 1 was 
more accurate than method 3. The pattern of the accuracy of 
the four methods was similar to that experienced in the 


Northern area. 
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fmeotrait of Juan De Fuca and Puget Sound 


Although method 4 demonstrated a significant increase 
in accuracy, none of the four methods which used historical 
salinity derived from the ICAPS file was deemed accurate 
enough. This file does not define any water mass applicable 
in the Strait of Juan De Fuca and Puget Sound, and the extrap- 
Blaeron of the ICAPS data to fit this region proved to be 
Motricient in accuracy. 

Peete: OF INFERRED SALINITY VALUE 

fee lH AND SIXTH METHODS) 

Salinity values of 31, 35 and 29 %o% were selected as 
representative of the Northern, Southern and Strait of Juan 
De Fuca/Puget Sound areas, respectively. These values were 
applied both throughout the entire water column (method 5) 
Bamecese were replaced by the observed surface salinity value 
tor the first layer only (method 6). The number of stations 
that exceeded the acceptable limits for each of these methods 
1s shown in Tables V and VI. A further analysis of these 
results follows. 

ine Northern Area 

wieene 1765 stations examined, 325 stations (18.33%) 
had a surface layer error that exceeded the acceptable limits 
when a value of 31 %» was used throughout the entire water 
monumm (method 5). Only 117 (6.6%) exceeded the limits in 


deeper layers, and only one in a layer deeper than 30 nm. 
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ABIES V. 


Pememo RESULTANT PROM THE USE OF A SINGLE VALUE OF SALINITY 


THROUGHOUT THE ENTIRE WATER COLUMN (METHOD 5)* 











Season: Wanter Spring Summer Fal} 

Avea Parircy ital S.b. <30 430 Total S.L. «50 290 TObah sol. <30 530 Toral S Ll. 5350 >30 

North Stevo 4 6cts COS y U 079 =—:1197 1u0 l 643 98 8 \. | 0 0 
(1307) 3.5) (0a) (29.0) (14.7) (0.1) (is .2 ele) (ee) (2-10) (0.0) (0,0) 

South S575 450 0 0 0 0 0 0 229.0 0 0 ly3 0 0 0 
(U.0) (0.0) (0.0) (U.0) (0.0) (0.0) (020) (0.0) (0,0) (0.0) (0.0) (0.0) 

Beeeresounds 299%, 127 22 16 3 149 38 16 15 100 928 Gz 18 Jo 25 15 i 

Strait of (1755) 9122-0) (6.5) C525) elOe 7)? (2061) (28.0) (2250) (18.0) (27570) °(1 5.5) (i295) 


Juiwr De Fuca 


*Table descriptions as in Table IV. 


TABLE VI 


ERRORS RESULTANT FROM THE USE OF A SINGLE VALUE OF SALINITY 


THROUGHOUT THE ENTIRE WATER COLUMN (METHOD 6)* 


— 























Season: Winter Spring a ponte es Fall 
Area Sulintry ae S.L. 530 >30 Total $.L. <30 230 Totul wi «tt >30 Total S... 350 >49 
ae Scena la Se oe asta eee 
" ears 22° an) a ee) 8 ore) “aug pens = fac a 0) We aa oni 
aa aay SP Oo oo 10.0) es om nil fueO) ut Keon t 0) iG) io} no) 
Priicor sae (0.0) (0.0) (4.7) ae (0.0) (0.0) a3; ee ey anal) as (0.0) (0.0) ts) 


Juun De Fuca 


a 





*Table descriptions as in Table IV 
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When the observed surface salinity value was used in the 
first layer, no stations exhibited an unacceptable error. 
Z. southern Area 
No excessive errors were found when either method 5 
Bro was applied, for which a constant value of 33 Vo was 
assumed for this area. 
Seeeotrait of Juan De Fuca and Puget Sound 
Pesalinitty value of 29 /66 was assumed for this area. 
Of the 472 stations examined, 113 stations (23.9%) exceeded 
the limits when this value was used by itself, but only 23 
feos) when it was substituted in the first layer by the 


observed surface salinity value. 


It has been seen that for the three areas, and for any 
Beethne methods considered, most of the errors exceeding the 
accuracy limits occurred in the surface layer only. This does 
not affect the accuracy standards when the depth is deeper 
than about 30 m, which represents most of the area over which 


hydrographic surveys are conducted. 
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Peco NCLUSIONS 


ijemresults shown in the previous chapter indicate that 
Metorical data can indeed be used in the computation of 
depth corrections due to variations in sound speed for the 
West Coast of the United States. Data derived from the ICAPS 
Salinity file are completely acceptable only in the Southern 
area, although good approximations can be made in the Northern 
area. The situation in the Northern area, however, can be an 
artifact created by the observed data format. If the surface 
value as observed in the sea is applied to adjust the histor- 
ical salinity curve from ICAPS, then the results are completely 
correct for all depths and seasons, both in the Northern and 
Bemenern areas. This by itself would indicate sufficient 
confidence to apply this method to actual hydrographic sur- 
veyS in the areas considered. It would only require the use 
of a table, such as Table I, and the measurement of the sur- 
face salinity, the temperature being measured with an XBT. 
The advantages of this method over what is presently being 
Seme are Significant, both in time and cost. 

ICAPS water masses are representative of very large areas 
Beeceep water. For application to coastal hydrography it is 
better to create a file of the shelf area exclusively. A 
ten year period may be representative, but a longer period, 


including as many observations as possible should be better. 
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Mme use of a single salinity value derived from a large 
number of observations 1s the easiest method to be applied 
mm the computation of depth corrections. One further step 
is to make a surface salinity observation and use that value 
for the computation of the first layer correction. The im- 
provement of this method over the method currently in use, 
from a practical point of view, is remarkable. 

mms Conclusion is applicable to the entire West Coast 
Seetcne United States, excluding the Strait of Juan De Fuca 
and Puget Sound where the results were not favorable, and 
to San Francisco Bay area where there were insufficient data. 

In summary, the applicability of the use of historical 
data in the computation of depth corrections is as follows: 

1. Northern Area 

A single value of salinity of 31 voeiseapplicable 

throughout the entire water column except for the surface, 
where the salinity must be measured. In these situations, 
the more complex method 4, wherein the historical salinity 
Meetile derived from the ICAPS' data file is adjusted to the 
surface value, may also be used. 

Z. southern Area 

The use of a single salinity value of 33 %,, appli- 

cable throughout the entire water column, is sufficiently 
accurate for this region. Alternatively, the historical 
Salinity curve derived from ICAPS can also be used with no 


adjustment for the surface layer. 


ao 





meottait of Juan De Fuca and Puget Sound 


tie Wise of hastonical Salinity data derived from any 


of the methods above 1S not recommended in this area. 
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APPENDIX A 


WILSON'S EQUATION 


Wilson first presented his empirical formula for the 
computation of the speed of sound in sea-water in the June 
1960 issue of the Journal of the Acoustical Society of Amer- 
ica [Wilson, 1960a]. In the October issue of the Same year, 
a Slightly modified version of this formula was introduced, 
Meaeieean extended range of applicability. This version is 
the formula commonly used and accepted as best representing 
the sound speed dependence on pressure, temperature and 
memrty. iit fits data within the ranges of -4°C to 30°C 
in temperature, l kg/em* to 1000 kg/cm? in pressure and 0 “jy 
memo?) /o in salinity, with a standard deviation of about 
oeou0 m/s {Wilson, 1960b] -. 

A value of 1449.14 is given for the sound speed at T=0°C. 
P=0.0 kg/cm? nc —95 5) pee OT recurons are added to this 
reference speed in order to compensate for the variations in 
meeoctre, temperature and salinity, either alone or in a 


combined effect. The revised formula is as follows: 


: An analysis of this formula for values of pressure, 
temperature and salinity outside the range for which the 
equation was derived showed that the error can go up to 
2.13 m/s for extreme conditions of temperature and salinity 
expected to occur in the sea [Wilson, 1962]. 
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V = 1449.14 + Ve a Vp + Vo + Vorp 


where V, = 4.5721T- 4.4532 x MONT omeoaas se ie” oT 
Pen ere 
= va arene 
V) = 1.60272 x 10> P + 1.0268 x 10°” P 
PGs 6h 5) S505 x 10 ~~ Po 
F aS 2 
V. = 1.39799(S- 35) + 1.69202 x 10°°(S- 35) 

; -2 -7 2 
Vorp* (S~ 35) (-1.1244 x 10° T+ 7.7711 x 10° T 
ae cient) Pee ods lo P 

ae aoe 
+ 3.1580 x 10 Phere. 5790 x 10 PT’) 
4 6 2 


ees oO7 ede ho 754812 x 10  T 


8 i 


Peer) - Foyer P-(22.5204 « 10” T 


10 


PMCS OSe Oe 91- eee (1 o646 x 10)” 7) 


The units of temperature, pressure, salinity and sound 
Speed are, respectively, degrees Celsius, kilograms per square 
centimeter, parts per thousand and meters per second. While 
temperature and salinity are usually directly measured, 
pressure 1s computed from depth. The dependence of pressure 
on depth is not linear and several expressions may be applied 
for this computation. An iterative process is used by NOS 
to determine pressure [{Wallace, 1971]. This method was 
extracted from the User's Guide to NODC's Data Services 
[NOAA, Department of Commerce, 1974], where it is described 


as follows: 
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somo sco sel + 8§.14876577 x 10° S 


- 4.8249614 x 10°" s* + 6.7678614 x 107° 5° 


B,980616 - 2.5928 x 10 “(cos29) + 6.9 x 10°° (cos29)° 


P(surface) = 10.1325 decibars 


g, + 1.101 x 10°" (d, +d,_,) 


ae 0 on)/R 
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A first approximation gives P’ and p;: 
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with Pl computed this way, p5 Poca! Gula eed. 


A second approximation using 05 gives Pr and O;: 
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with Pr 5 is calculated. 
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fee finally, with ep; 
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Mmimene above expressions, the subscript indicates the depth 
mudeeis being computed, and the number of primes is inversely 
[eoportional to the accuracy of the approximation. The mean- 


ing of the different parameters is the following: 


d 


depth in meters 

T = temperature in degrees Celsius 

Ew Salinity in parts per thousand 

o = density 

fae hessunre in kilograms per square Centimeter 


e- slatitude in degrees. 
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Figure B-1. Location of stations for the northern area 
and Strait of Juan De Fuca/Puget Sound, in winter. 
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meeune B-2. Location of stations for the northern area 
and Strait of Juan De Fuca/Puget Sound, in spring. 
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Figure B-3. Location of stations for the northern area 
and Strait of Juan De Fuca/Puget Sound area, in summer. 
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Figure B-4. Location of stations for northern area and 
strait of Juan De Fuca/Puget Sound area, in fall. 
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Figure B-5. Location of stations for southern area, winter. 
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meeure B-6. Location of stations for southern area, spring. 
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Figure B-7., Location of stations for southern area, summer. 
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Figure B-8. Location of stations for southern area, fall. 
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SUBROUTINE STAT (M,AsMEAN,STODEV) 


SUBROUTINE STAT COMPUTES THE MEAN AND STANDARD 
DEVIATION OF THE ARRAY Ay WHICH HAS M ELEMENTS 


DIMENSION MUST BE ASSIGNED IN MAIN PROGRAM FOR ARRAY A 


DIMENSION A(M) 
REAL MEAN 


COMPUTE MEAN 
IF nen GTe-1) GO TQ 800 
MEAN=A (1) 
STDEV=0. 
RETURN 
800 SUM=0. 
DO 801 [=1,M 


SUM=SUMtA(T) 
801 CONTINUE 
ME AN=SUM/FLOAT(M) 


COMPUTE STANDARD DEVIATION = 


) -MEAN ) *¥*2 
ERM 


802 
UM/FLOAT(M=-1)) 
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